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Abstract: We compute the two photon exchange contributions to elastic scattering of polarized elec- 
trons from target protons. We use a nonlocal field theory formalism for this calculation. The formalism 
maintains gauge invariance and provides a systematic procedure for making this calculation. The results 
depend on one unknown parameter b. We compute the two photon exchange correction to the ratio of 
electric to magnetic form factors extracted using the polarization transfer experiments. The correction 
is found to be small if & 1. However for larger values of 6 > 3, the correction can be quite significant. 
The correction to the polarization transfer results goes in the right direction to explain their difference 
with the ratio measured by Rosenbluth separation method. We find that the difference between the two 
experimental results can be explained for a wide range of values of the parameter b. We also find that 
the corrections due to two photon exchange depend on the photon longitudinal polarization e. Hence 
we predict an e dependence of the form factor ratio extracted using the polarization transfer technique. 
Finally we obtain a limit on b by requiring that the non-linearity in e dependence of the unpolarized 
reduced cross section is within experimental errors. 
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1 Introduction 



The observed discrepancy [2 [2] in the proton electromag- 
netic form factors in the JLAB polarization transfer ex- 
periments [31 m El E] and the SL AC Rosenbluth separation 
experiments [3 [HI E] has been studied extensively in the lit- 
erature. The two photon exchange contributions are most 
likely the source of the difference [TUl El [H [H [H [H] . 
In a recent paper we constructed a nonlocal Lagrangian 
to model the electromagnetic interaction of extended ob- 
jects such as the proton pjjj. The model maintains gauge 
invariance in the presence of a form factor at the electro- 
magnetic vertex. We truncate the Lagrangian to include 
only operators with dimension five or less. The dimension 
five operator is necessary if we include the contribution 
proportional to the Pauli form factor F2. This truncation 
is reliable as long as the off-shellness of the proton propa- 
gator is small compared to the hadronic momentum scale. 
The resulting Lagrangian depends on the two on-shcU form 
factors, Fi and -Pj, and contains one unknown parameter 
which we denote as b. Wc found that for small values of this 
parameter 6 1, the results of the SLAG Rosenbluth mea- 
surement, after the two photon exchange correction, are in 
agreement with the JLAB polarization transfer results. In 
the present paper we compute the two photon exchange 
corrections to the polarization transfer experiment. 

In polarization transfer experiments, longitudinally po- 
larized electrons are scattered from fixed target protons. 
The tree level amplitude for the process (See Fig. [T]) is 



given by: 
where 



--u{k')-|^'u{k, Se)U{p', Sp,)r^U{p), (1) 



Fiiq)l, + -^F2iq)a,^q". 



2Af, 



(2) 



Here e is the charge of proton, fc, k' are the momenta of 
the initial and final electron, p, p' are the initial and final 
proton momenta, Mp is the proton mass, Kp = 1.79 is pro- 
ton anomalous magnetic moment and q = p'— p = k — k' 
is the momentum transfer. 

Using ([T]), one can calculate the tree level cross-section 
for the scattering of longitudinally polarized electrons from 
protons: 



where |A^oP represents the amplitude squared, summed 
over final electron spin and averaged over initial proton 
spin and is given by: 
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^ tr[{(^ + m,) (1 + -751^)7'' + ^n^h') 
+ A/p)7orj:7o(/ + Mp) (1 + 75^pOr, 



X tr 



(4) 
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polarized recoiling proton. Let us define, 




P 



Figure 1: The one photon exchange diagram contributing to 
the elastic electron proton scattering. Here k, k' refer to the 
initial and final electron momenta and p, p' to the initial and fi- 
nal proton momenta respectively. The symbol q = k — k' = p'—p 
denotes the momentum exchanged. 

Here s^, is the spin four-vector of the final proton. In 
general the spin four-vector s for a particle with mass m 
and momentum p = {E; p) can be written in terms of the 
unit three- vector n, specifying the spin direction in the rest 
frame of the particle, by 



n ■ p ^ n ■ p 
; n + p- 



m{m + E) 



(5) 



For the incident electron, energy E^, is much greater than 
its mass me, and we have approximated the spin four vec- 



tor, w hkf^/me, where h 



■ k. 



Let the scattering plane be the X — Z plane where the 
momentum of the recoiled proton p' defines the Z axis, 
i.e., z = p' . The Y axis is defined a,s y = k x k' . This 
also defines the X axis a.s x = y x z. With this choice of 
coordinate axes, the tree level cross section can be written 
as: 



where 




{hp, ■ z) + hl^'^ P^'^ {hp, ■ x) , (6) 



jl7 _ 



rl7 pl7 _ 



7-17 pl7 



Ee + E' 



v/r(l+T)G^tan2^, 



-2VT(l + r)GMGEtan^. 



Here r = -q^/AM^, e = 1/[1 + 2(1 + r) tan2(6le/2)] is the 
longitudinal polarization of the photon and 9e is the elec- 
tron scattering angle, /q'^ is proportional to the tree level 
unpolarized cross-section. Iq'^ P^'^ [I^'^P^'^) is proportional 
to the term in the tree level cross-section that corresponds 
to scattering a longitudinally polarized electron from an un- 
polarized proton producing a longitudinally (transversely) 



d{Aa}p) 
dQp 



a'^E'^ cos2 ^ 



,4 



^ 8Ei sin^ ^ 

' a'^E'^ cos^ ^■ 
8EI sin^ % 
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1 + r' 

7I7 pl7 
-'O 

1+T 



7-I7 pl7 

1q T-rj, 



(7) 
(8) 
(9) 



Using these we can rewrite da^'^ / dQ,e as 



d<7, 



17 



dQ.e 

d{/^a}g 



d{Aal'') 



+ 



drip 



dile 

'{hp' ■ x) 



{hp, ■ z) 



(10) 



Hence, the ratio of the form factors is given by: 



R 



Ge 
Gm 



P^^ Ep + E' 



pl7 



2M„ 



■ tan — 
2 



d{Aa^^)/dnpEe + E '^ e. 
'd{Aa]^)/dnp 2Mp 



tan — . 
2 



(11) 



The contribution of the two photon exchange diagrams 
to the electron-proton elastic scattering cross section can 
be written as 



dcr27 ^ 2Re{M*oM2'y)E, 



'2 



dQp 



GAM^TT^E^ 



+ 0{a'), 



(12) 



where A^27 is the total amplitude of the two photon ex- 
change diagrams. As in the case of tree level process we 
can define d{Aa'^)/dne and d{Aa'^)/dfle from the terms 
in da'^'^ /dftp that are proportional to h{hp' -z) and h{hp, -x) 
respectively. Then, 



da^ 
dClp 



da. 



dn 

+h 



dQp 



■{fip' ■ y) 



d{Aa2_ 
d{Aa^) d{Aa?g) 



dflp 



-{hp, ■ z) + 



dflp 



{hp, ■ x) . (13) 



Here da^^^ /dflp represents the terms independent of the 
spin of the final proton and d{Aa'^)/d^le corresponds to 
the normal polarization of the final proton. 

Experimentally the polarization components, Pl and 
Pt arc measured for different q^. From these the ratio of 
the form factors is extracted by the following relation: 



Pt Ep+E' Op 

tan — 

Pl 2Mp 2 



(14) 



In the tree level approximation Pe is identical to R. Let 
Gm be the experimentally measured value of Gm- Then, 
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(a) 



(b) 



(c) 



Figure 2: The two-photon exchange diagrams contributing to the elastic electron proton scattering: (a) box diagram, (b) cross- 
box diagram and (c) diagram proportional to . Another diagram proportional to P is obtained by interchanging / <-» (g — Z) in 
(c) 



with the knowledge of Gm and Re one can define the fol- 
lowing quantities: 



d(A(TT) 



a^E' cos^ 



^El sin^ ^■ 



' a^E'cos^^ 



8E^ sin' 



4 



1 + r' 

IqPt 
1 + r 



where 



IoPl 



IoPt 



E,+E' 



-^y7(TT7)GLtan^^, 



t)GIiRe tan- 



(15) 
(16) 

(17) 
(18) 



and (b) show the box and cross-box diagrams respectively 
and Fig. [^c) shows a diagram proportional to Ir. The 
amplitudes for the box and cross-box diagrams are given 
by: 



iMb = 



From these one obtains the ratio R corrected for the two- 
photon exchange processes by the following relation: 



R 



d{AaL)/dne - d{Aal^)/dne 



E, 



Re 



tan — 

2A/p 2 

1 - At 



1- At 



where 



d{Aa^)/dVl 



At 
Al 



d{AcTT:)/dne 



d{Aa^)/dVle 



d{Aai^)/dne 



(19) 
(20) 

(21) 
(22) 
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2 Calculations and Results 



dH 



e(fcV(/^-/)7'-M(fc,Sfc) 
{{k - l)^ ^ ml + i^) 

1 



U{p\ v)<i FMl^^ + ^^F2{q)a^^r 



' + / + Mp 



X i^l(07. + *^i^2(0^./3^^ C/(p) 



, (23) 



dH 



[k-iy- ^ml + iS, 
1 



2Af„ 



(24) 



As explained in |16j there arc three diagrams that con- 
tribute in the two-photon exchange processes. Figs. [DJa) 



Here q = q — I and ^ is an infinitesimal positive parameter. 
The amplitude for the diagram proportional to can be 
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Figure 3: Contribution of box and cross-box diagrams to Al (unfilled circles) and At (filled circles) for different - (a) 
= -2.47G•eF^ (b) = -2.97Ge"l/^ (c) = -3.50Ge\/^ (d) q^ = -3.97Ge\/^ (e) q^ = -4.75GeV^ (f) q^ = -5.54GeF2. 
At is defined in (l2T|l and Al is defined in (l22ll . 
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Figure 4: Contribution of the diagrams proportional to to Al (unfilled circles) and At (filled circles) for different - (a) 
= -2A7GeV^, (b) = -2.97GeV^, (c) q^ = -S.50GeV^, (d) = -S.97GeV^, (e) q^ = -4.75GeV^ (f) g^ = -5.54GeV^. 
At is defined in (l2T|l and Al is defined in (l22ll . Here b = 1. 
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written as: 



8M2 



dH 



iZ(fc')7''(^-/)7Mfc,Se) 



The total two photon exchange amphtude, 



M 



27 



Mb + Mcb + Mj, 



(25) 



(26) 



In the numerical calculation we have dropped mg. The 
model used for the form factors [17l [18] is the same as 
Model I described in (TB]. The values of the parameters 
are repeated here in Appendix A for convenience. The re- 
sults of our calculation may show some dependence on the 
precise form factor model. This appears to be the main 
uncertainty in our calculation which we hope to explore in 
a future publication. 

Contributions from box and cross-box diagrams are 
computed at 10 different values of /i^ (from 0.005 to 0.0095) 
for each value of and e. The cross-box diagram is well 
defined but for the evaluation of the box diagram a small 
imaginary term ^ is kept in the propagators. This makes 
the integral in the infrared limit well defined in the case 
me = 0. For each value of g^, e and jj.^ we have calcu- 
lated the box diagram amplitude for 4 different values of 
^ (between 0.001 and 0.00175). The final dependent 
box diagram amplitudes are obtained by extrapolation to 
^ = 0. The diagram proportional 6^ has no infrared (IR) 
divergent term. So the contribution coming from it is com- 
puted keeping = 0. 

Following Mo and Tsai [HI EQI Ell [22] one can show 
that in the limit —> the leading term from the box and 
cross-box diagram can be expressed as: 



(27) 



where 



dx 



{^Pi + (1 - x)pjf 



X In 



[xpi + (1 - x)pjf 
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Figure 5: The ratio UpGe/Gm from polarization transfer ex- 
periments (filled circles) and after correction for the two-photon 
exchange diagrams with 5 = (unfilled squares) ,b = 2 (unfilled 
circles), b = 3 (unfilled triangles), b = 4 (unfilled diamonds), 
b — 5 (unfilled pentagon) and b = 6.91 (filled triangles). The 
corrected values have been offset by AQ^ =0.1 GeV^ for clar- 
ity. The filled diamonds and the filled squares represent the 
Rosenbluth extraction experiment results obtained at JLAB and 
SLAC respectively. 
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Figure 6: The ratio fipGE/Gm from SLAC Rosenbluth extrac- 
tion experiments (filled squares) and after correction for the 
(unpolarized) two-photon exchange diagrams with b = (un- 
filled squares), 6 = 2 (unfilled circles), 6 = 3 (unfilled trian- 
gles), 6 = 4 (unfilled diamonds), 6 = 5 (unfilled pentagon) and 
6 = 6.91 (filled triangles). The filled diamonds represent the 
JLAB Rosenbluth extraction experiment results and the filled 
circles represent the results of the polarization transfer experi- 
ments. 

Hence, the leading contributions to the cross sections com- 
ing from the box and cross-box diagrams are given by: 



dflp 



2a 



K{p',k)-K{p,k)] 



IR 



d(A<y 

dfle 



(28) 
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Figure 7: Dependance of the form factor ratio ^pGe/Gm on e from different - (a) = -2.47GeV^, (b) q'^ = — 2.97GeV^, (c) 
q^ = -3.50GeF^ (d) = -3.97GeV^, (e) q^ = -4.75Ge'l/^ (f) = -5.5AGeV^. For each q^, the ratio, fiG^/GM is plotted 
for different b values -6 = (unfilled squares), 6 = 2 (unfilled circles), 6 = 3 (unfilled triangles), 6 = 4 (unfilled diamonds), 6 = 5 
(unfilled pentagon) and 6 = 6.9f (filled triangles). 
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Let 



-'L,T 



In ^2 



(29) 



IR 



To remove the IR part from d{Aa'^rp)/dne wc fit these 
with the following funetions: 



with 



d(Aagr) 



&l,t(a^') 



aL,T(M^) + foL,T(M^)lnM^ 



(30) 



12 



L,T M 



Hence, x give the IR removed d{A(7-^rj.)/dCle- 

Fig. [3] shows the contribution to Al,t coming from box 
and cross-box diagrams for different q^. The contribution 
from the diagram proportional to is shown in Fig. U) 
Fig. [Hlshows the ratio fipGE/Gu from polarization transfer 
experiments (filled circles) and after correction for the two- 
photon exchange diagrams with 6 = (unfilled squares), 
b ~ 2 (unfilled circles), 6 = 3 (unfilled triangles), 6 = 4 
(unfilled diamonds), 6 = 5 (unfilled pentagon) and 6 = 6.91 
(filled triangles). The highest value of 6 used is obtained 
from the limit on nonlinearity of JLAB Rosenbluth data 
(see Section 3). In the range 3.5 > > 2.5 GeV^, the po- 
larization transfer experiment used the kinematic regime 
with e lying roughly between 0.7 and 0.8 |5j- For simplic- 
ity, here we assume a fixed value e = 0.75. Fig. [5] shows 
the effect of 6 on the unpolarized SLAG data. The method 
used to obtain these corrections is described in Section 6 of 

The form factor ratio Ge/Gm depends only on Q^. 
However the ratio i?E, defined in (jl4p . may show a de- 
pendence on £ due to higher order contributions. We pre- 
dict this dependence by assuming that the electron-proton 
elastic scattering cross section can be well approximated 
by including only one and two photon exchange diagrams. 
The uncorrected form factor ratio i?E can be calculated by 
using 1^0]) . where the corrected ratio R is independent of 
e. Fig. [7] shows the predicted dependance of the uncor- 
rected form factor ratio ^pi?E on e for different 6 and q^. 
Here again we have assumed that the kinematic regime of 
the polarization transfer experiment corresponds to a fixed 
e — 0.75. This may be used in future to fix the value of the 
parameter 6. 



unpolarized electrons from protons, ctr, with the following 
function: 



(31) 



where e = e - i. As the resuh obtained by JLAB [251121] 
contains much smaller error-bars than the SLAG data we 
use JLAB data to obtain V2 for different = —q^- Then 
we fit the contribution of the box and cross-box diagrams 
(cr^'-'^) and the contribution proportional to 6^ (cr^) to the 
reduced cross-section with similar functions, i.e., 



,BCB 



b^{nl+TZle 



BOB; 



(32) 
(33) 



To put the limit on 6 we first obtain a limit on the 
nonlinearity parameter 7^2- The 2 sigma limit is obtained 
by finding the largest value of —V2 such that deviates 
from its minimum value by 4 units. We consider the largest 
value of —V2 since the two photon exchange contributions 
are found to give a negative value of 'p2- The parameters 
Vo, Vi are also allowed to vary while determining the max- 
imum value allowed for —V2- The best fit value of V2 and 
the 1 and 2 sigma limits are given in Table [T] The best 
fit values of TZ^'^^ and are given in Table [H Assuming 
that the dominant source of nonlinear behaviour are the 
box, cross-box and 6 diagrams, we obtain the limit on Ir 
by the following relation. 



'P2limit ~ 7?.f 

7^1 



(34) 



where V2Umit is the one or two sigma limit on the parameter 
7^2- From Tabled] we see that the most stringent limit on 6 
is obtained for = 2.64 GeV^ and is given by: |6| < 6.91. 
The nonlinearity of the reduced cross section ctr has also 
been computed in |25] using the GPD formalism. 



r2 V2 V2 62 

in best 1 sigma 2 sigma 2 sigma 
GeV^ fit limit limit limit 



2.64 2.02 -3.82 -9.64 47.74 
3.20 1.10 -3.70 -8.52 70.61 



3 Limit on h 

As described in [16] 6 is a free parameter. But the fact 
that (T^ varies nonlinearly with e allows us to put a limit 
on 6. We fit the reduced cross-section for the scattering of 



4.10 2.68 -2.92 -8.48 152.93 



Table 1: The best fit and the one and two sigma limit on the 
parameter 7^2 (defined in (|31[) ) for different . All values of V2 
have been scaled by 10'*. The corresponding two sigma limit on 
the parameter I? is also given. 
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Q2 in GcV2 



A.. 



2.64 -4.04 X 10-"* -11.73 x IQ-^ 
3.20 -3.69 X lO-"* -6.84 x 10"^ 
4.10 -1.95x10-* -4.27x10-6 



Tabic 2: Values of the fitting parameters T^.f'^^ (defined in 
|)) and TZl (defined in ^) for different . 



4 Conclusion 

We have computed the two photon exchange corrections 
to the proton electromagnetic form factor ratio ^ipG^/ Gm 
using a gauge invariant nonlocal Lagrangian. The La- 
grangian is truncated to dimension five operators and de- 
pends on one unknown parameter h. The higlier dimension 
operators are expected to give negligible contributions as 
long as the ofF-shellness of the proton propagator is small. 
The two photon exchange corrections to the ratio are found 
to be small as long as the parameter b ^ \. We impose a 
limit on this parameter by tlie predicted nonlinearity in the 
e dependence of the unpolarized reduced cross section ctr 
due to two photon exchange contributions. A two sigma 
limit is found to be |&| < 6.91. For such large value of h 
the corrections to the polarization transfer experiment are 
quite significant. The corrections are found to go in the 
right direction and make the polarization transfer results 
come very close to the SLAG Rosenbluth measurement of 
this ratio. Hence we find that for a wide range of values of b 
the results of the two experiments agree within errors after 
including the two photon exchange contributions. We also 
predict an e dependence of the ratio extracted from polar- 
ization transfer. This can be tested in future experiments 
and also be used to extract the value of the parameter h. 



Appendix A Model for the Form Factors 

The fits for GM/^J'p and Ge are given by: 
Guiq^) A A 



fj-p 



a=l 



GE(g^) = E 



^ (g2 -ml + imaV'J ' 



(A.l) 
(A.2) 



The values of the masses and the parameters are tabulated 
in Table[3l Using the fits for the magnetic and electric form 
factors one can determine the Dirac and Pauli form factors. 



1 


-2.882564 


-3.177877 


0.8084 


0.2226 




-hi 1.944314 


-fi 2. 123389 






2 


2.882564 


3.177877 


0.9116 


0.1974 




-i 1.944314 


-i 2.123389 






3 


-1.064011 


-0.608148 


1.274 


0.5712 




-^3.216318 


5.685885 






4 


1.064011 


0.608148 


1.326 


0.5488 




-fi 3.216318 


+i 5.685885 






5 





3.211388 


1.96 


1.02 






-Fi 0.693412 






6 





-3.211388 


2.04 


0.98 






-i 0.693412 







Table 3: Masses, widtlis and parameter values for Gu/l-ip and 
Ge fits. A s and B s are defined in (|A.ip and (|A.2p . 
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